Order parameter statistics in the critical quantum Ising chain 
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In quantum spin systems obeying hyperscaling, the probability distribution ol the total magne- 
tization takes on a universal scaling form at criticality. We obtain this scaling function exactly for 
the ground state and first excited state of the critical quantum Ising spin chain. This is achieved 
through a remarkable relation to the partition function of the anisotropic Kondo problem, which 
can be computed by exploiting the integrability of the system. 



The concept of universality is central to our under- 
standing of continuous phase transitions. Universal phys- 
ical quantities at or near a transition in one system co- 
incide with those in other systems that share a few key 
characteristics, including dimensionality and symmetry 
of the order parameter. Continuous transitions are there- 
fore grouped naturally into universality classes of com- 
mon critical behavior pQ. Familiar universal quantities 
include critical exponents that characterize the singular 
behavior of thermodynamic and response functions in the 
vicinity of the critical point . Although the amplitudes of 
these singularities are entirely system-dependent, certain 
amplitude ratios are universal [2]. 

A particularly natural family of universal amplitude 
ratios is formed from the volume-integrated order pa- 
rameter M in a finite-size system by 

a {M2n) m 

It is readily seen that the hypothesis of universality ap- 
plied to this entire family is equivalent to the existence 
of a universal scaling function f(X) defined by 

f(X) = sP(sX) (2) 

where s 2 = (M 2 ) is the variance of the order parameter 
and P(m) its probability distribution function. The exis- 
tence of this scaling function relating the probability dis- 
tribution for different system sizes is a consequence of hy- 
perscaling, as emphasized in Ref. 3. Extensive numerical 
work has confirmed the universality of sP(sX) [4j[5]. For 
the important benchmark of the two-dimensional classi- 
cal Ising model, Ref. [HI is a useful guide to the literature 
as well as a tour-de-force numerical study. 

The aim of this work is to remedy two striking defi- 
ciencies in our present understanding of order parameter 
distributions: the small number of analytical results and 
the scant attention that has been paid to quantum phase 
transitions. The study of order parameter fluctuations 
at quantum critical points, aside from being inherently 
interesting, is motivated by recent experiments in atomic 
physics, where the measurement of the full distribution of 
global fluctuating observables has become a reality [7J [H] . 



For the quantum Ising chain, we compute P(m) exactly 
by relating its generating function to the partition func- 
tion of a particular anisotropic Kondo problem. This re- 
markable relationship allows the application of the pow- 
erful analytic methods developed to solve this and other 
quantum impurity problems. 

The quantum Ising chain, often referred to as the 
transverse-field Ising model, has Hamiltonian 

L 

H=-Y,[h<Tf + o-!*! +1 ] . (3) 
i=l 

For the moment, we impose periodic boundary conditions 
so that cr£ +1 = a\. The critical point of the model Eq. (l3| 
is at h = 1, separating the ordered (h < 1) and disordered 
(h > 1) phases. The magnetization M = J2i a i/^ does 
not commute with H, so eigenstates of H are typically 
sums over states with different eigenvalues of M. The 
distribution functions of the magnetization in the ground 
states for various boundary conditions were studied nu- 
merically in Ref. HJ and good scaling of sP(sX) onto a 
universal curve was found for L > 16. There are also 
analytical results for the distribution of the transverse 
magnetization of [51 1X0]. 

The magnetization distribution in a state |a) is 

P a (m) = (a\8(m- M)\a). 

Rewriting the delta function as an integral gives 

/°° rl\ 

where x Q (A) = (a|e lAM |a) is the generating function 
of the moments of the distribution. The flip operator 
T = Yii o~f commutes with the Hamiltonian, so the re- 
sulting Z2 symmetry requires that P a {m) = P a (—m). A 
generalized Lee- Yang theorem shows that the generating 
function for the ground state has the factorization [IT] 

*„ ( A)=n(i-f) (4, 
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for real positive {E p }. The even cumulants, given by the 
coefficients of the expansion of In xo (A) , are then 



<M 2 ") C = (-1) 



,-i(2n)! 



and therefore alternate in sign. 

In the scaling limit, the sum in the operator M can be 
replaced with an integral M = dxcr(x), where a(x) is 
the standard Ising quantum field. Expectation values in 
the ground state can be computed in the path-integral 
picture by taking Euclidean spacetime to be a very long 
cylinder of circumference L\ the long cylinder means that 
in the Hamiltonian picture the system is projected onto 
its ground state. The path integral for the generating 
function for the ground state is then precisely the par- 
tition function of the 2d classical Ising model with an 
imaginary magnetic field along a defect line wrapping 
around the cylinder. In this path integral, we are free 
to exchange the roles of space and time so that the new 
Euclidean "time" direction r is periodic. Since all the 
operators in xo were originally at the same time, this ex- 
change puts them all at the same spatial position. Thus 
in this new picture, xo describes the continuum limit 
of infinitely-long Ising chain at temperature 1 / L and an 
imaginary magnetic field iX at point 0. The underlying 
lattice Hamiltonian is that of the quantum Ising chain ^ 
on an infinite line, so that the sum runs from i — — oo to 
oo. The generating function is simply 



XoW =tr 
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(5) 



where T represents time-ordering and the (Euclidean) 
time dependence of the operator denotes the Heisenberg 
picture: (Jq{t) = e Hr <jQe~ HT . 

Writing Xo as (J^J allows us to compute it exactly at the 
critical point. We first relate Xo to the partition function 
of a famous quantum impurity model, the anisotropic 
Kondo problem [12] . We then apply the methods of in- 
tegrability to compute XoW an d hence P^m). 

There are two ways of showing why the Kondo prob- 
lem arises. The first is quite direct. At the critical point 
h = 1, the trace in Eq. §5§ can be expressed as an ex- 
pectation value in the critical Ising field theory. We then 
can expand in A and use the known spin correlation func- 
tions in the field theory to write integral expressions for 



the moments. The integrals for the term order A 2 ™ term 
are over values of < Tj < L, but because spin correla- 
tions are independent of the ordering of the Tj, we can 
multiply by (2n)\ and order them < t\ < ■ ■ ■ r 2 „ < L. 
This allows us to exploit a marvelous result for critical 
Ising correlators when all spin fields lie on a cycle of a 
cylinder [T3] : 



(d(0,T 2n )---<r(0,Ti)) oc J] 



i> 3 
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(6) 

where g — 1/8, the scaling dimension of the spin field. 
This formula is valid only when the Tj are ordered. 
Keeping track of the constants in front of ^ shows 
Xo(A) is a function of the dimcnsionlcss quantity z = 
A(27ra) 1 / 8 L 7 / 8 , which establishes the system size inde- 
pendence of the amplitude ratios Eq. ([lj and the scaling 
form of the order parameter distribution (Eq. ^ ) for this 
model, with s cx a 1 / 8 ^ 7 / 8 . 

The correspondence with the Kondo problem is now 
apparent. Xo(A) takes the form of the partition function 
at imaginary fugacity iz of a Coulomb gas on a ring. The 
gas consists of positive and negative charges with loga- 
rithmic interactions; because of the (— l) t+J in ^ the 
signs of the charges required to alternate in space. This 
is precisely the celebrated Anderson- Yuval expansion for 
the partition function Zk of the Kondo model, describ- 
ing the interaction of a spin- 1/2 impurity with a Fermi 
gas |12j . The expansion in z is a perturbative expansion 
in the spin- flip part of the Hamiltonian, with alternation 
in signs of charges arising from the two spin states of 
the impurity. In the Kondo problem g parametrizes the 
anisotropy (in spin space) of the interaction. Thus we 
have shown that 



Xo(A) = -Z K (iz), 



(7) 



with g = 1/8. A very similar result was derived for the 
chiral version of this problem, describing a point contact 
in a p + ip superconductor |14j , and is in accord with the 
results for the boundary entropy of the Ising model with 
a defect magnetic field [151116]. 

An illuminating way of rederiving Eq. ^ is to use 
boundary conformal field theory. Expanding Eq. ^ in 
A and using the time ordering gives 
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(8) 



Since a z a x a z — —<r x , the effect of the boundary mag- netic field is to flip the sign of the transverse field h at 



3 



1.5- 

Z K (i\) : 
1- 

0.5- 



-0.5 



-n 



FIG. 1: Zpc(iX) at g=l/8 and (dashed) its approximation (101 




site zero for every other interval between insertions of 
<7g. This means we can effectively treat h as being time- 
dependent, i.e. h(r) = h JljZi s E n ( T j ~ T )- The Ising 
spin chain has only nearest-neighbor interactions, so we 
can "fold" the theory in half at site 0, turning the de- 
fect into a boundary. The continuum analog of CTq (t) is 
then a boundary- condition- changing operator B{t) |17j . 
Here, inserting B(tj) at an instant Tj toggles between two 
different critical boundary conditions. The bosonization 
analysis of Ref. [15] allows us to show that B(t) is iden- 
tical to the spin-flip operator in the anisotropic Kondo 
model. This follows from two key facts about B(t): it has 
dimension 1/8, and inserting it simply toggles back and 
forth between two Dirichlet-type boundary conditions on 
the boson, corresponding to fixed values ip = w/4 and 
ip = 3tt/4 of the boson field. A free boson calculation 
immediately yields Eq. ^ with 2g — A</? 2 /7r 2 = 1/4. 

We now find XoW> its asymptotics, and its moments. 
There are three distinct ways of evaluating the partition 
function Zx{iz), all of which work to high numerical ac- 
curacy. One way is to use the thermodynamic Bcthc 
ansatz [18] H5] . the second is to use series expansions 
[HI [2Q], and the third is to compute the spectral de- 
terminant of an associated ordinary differential equation 
[5TJ[55]. We have used the latter two approaches, both of 
course giving the same result, displayed in Fig. [I] along 
with a very accurate asymptotic expression to be dis- 
cussed below. In accord with the generalized Lee- Yang 
theorem mentioned earlier, the zeroes Z^{iz) occur at 
real z. In fact, the set {E p } appearing in Eq. Q are just 
the eigenvalues of the spectral problem. 

The scaled distribution function sP(sX) is obtained 
by (numerically) taking the Fourier transform of Zx(iz), 
and shown in Fig. [2] We have tested our predictions by 
exact diagonalization of the lattice Hamiltonian Eq. ([3| 
for L = 16 — 23 using the ALPS libraries [23] HJ. We find 
excellent agreement, with the numerical results getting 
closer to the exact curve as L is increased. 



FIG. 2: Dots show numerical data for sP(sX) for the ground 
state of the critical Ising chain for L = 16 — 23. The red line 
is the scaling function fo(X) for the ground state, the Fourier 
transform of Zn(iz). 



We find accurate approximations for the scaling func- 
tions and exact values of the moments by utilizing the 
functional relation between the partition functions of the 
Kondo and the boundary sine-Gordon models [TS] [55] 156"] 
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-Zbsg(z) 

Both partition functions have similar Coulomb-gas ex- 
pansions, but Zbsg is given in terms of unordered 
charges. The asymptotic expression for Zbsg{z), valid 
for large z in a region near the positive real axis is |19j 

1/(1-9) 

i / 1 H III / illl \ 
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where a = ng/(2 — 2g). (In the differential equation 
approach [5TJ [55] , this expression arises from the WKB 
approximation.) Plugging this into (|9| gives 

-sin(a) (rnoA) 1 /' 1 " 3 ' 



Z K (iX) ps 2 cos cos(a) (m A) 1 -<> 



(10) 



where 



to 



2 v^ r (i-g) \ r(i-,g) 
r(i-f) J ^ 



As A — > oo along the real axis, this falls off as an oscillat- 
ing exponential. Even though this asymptotic expression 
in principle holds only at large A, when g — 1/8 it is a 
quite accurate approximation. It is plotted along with 
the exact expression in Fig. [T] Taking the Fourier trans- 
form of the asymptotic expression, one finds 
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FIG. 3: Comparison of numerical and analytical results for 
the first excited state scaling function fi(X). 



while a stationary phase approximation to the Fourier in- 
tegral gives fo(X — > oo) oc 

x - i+i/2 g ex p(-cX 1 /s) with 

c=g[m /(l-g)} 1 - 1 ^. 

The moments of the distribution are related to Zk by 

(M 2n ) oc (2n)\K 2n /2, where 

OO OO 

Z K {z) = K 2n z 2n , Z B sg(z) = B2nz 2n . 

ri— 1 n—1 

The universal amplitude ratios are therefore 

(2n)! K 2n 
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The B 2n have explicit series expressions [5D] which allow 
accurate numerical evaluation; the K 2n are then found 
in terms of the B 2m with m < n using Eq. ([9]). The 
first two ratios are exactly A4 = 1.43138 . . . and As — 
2.35464 . . . , while from exact diagonalization we find 
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This method can be extended to the computation of 
the distribution function for excited states as well. The 
lowest-lying excited state can be viewed as the ground 
state of a periodic Ising system with an antiferromag- 
netic defect. In the Hamiltonian formulation of the field 
theory, this excited state is given by acting on the ground 
state with the spin field, so the generating function for 
the moments of the first excited state is given by insert- 
ing spin fields a(x — +00, r) and a(x = — 00, r) into 
the correlator These insertions amount to setting 

p = 1/4 in the conventions of [20 or p — 1/8 in [25] . 



The resulting scaling function fi(X) of the first excited 
state is shown in Fig. [3| It is quite different from that of 
the ground state, illustrating the importance of different 
boundary conditions. 
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